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HARMONIC ANALYSIS OF DISTORTION IN AMPLITUDE MODULATION SYSTEMS 

WITH ASYMMETRICAL SIDEBANDS 



SUMMARY 

The general case of two modulating tones is mainly considered, and the 
amplitudes of the various frequency components appearing at the output of an 
ideal envelope detector are evaluated, including cases with unequal upper and 
lower sidebands. Difficulties associated with overmodulation are discussed, 
and a general method of finding the magnitude of any distortion product is ex- 
plained. This method is not affected in a particular case by the modulation 
depth. 

Coefficients associated with the lowers-order distortion products have been 
evaluated (Table 1) for cases in which the sidebands have the ratio 1, V2, 1/3 
or 1/9, and for the single- sideband case. For comparison between theoretical 
and experimental results, a measurable quantity has been calculated; agreement 
is remarkably good. 



1. INTRODUCTION 

Ideally, an analysis of a multltone modulating 
system of n tones would be valuable, but the com- 
plexity of such an analysis increases violently with 
n. The case n = 1 is special, and in what follows we 
have therefore taken n to be 2, and given special 
atte'ntion to cases where there is symmetry, or any 
other feature which simplifies the analysis. 

The main difficulty of the analysis is that the 
output from an envelope detector is such that only 
its magnitude, S{/), is determinate; the detector out- 
put is essentially positive* It is necessary to take a 
square root in the process of determining S{t)t and the 
positive square root must be chosen. Now the positive 
square root of x^ is not x unless we are certain that x 
is always positive; in general, the positive square 
root is lx|, and this expression has an awkward sing- 
ularity when X is capable of changing sign. In the 
case of S{t), this difficulty can arise whenever there 
is over-modulation. 

In Section 2, S(t) Is formulated implicitly for the 
general case, while in Section 3 certain cases when 
S{t) is simple and explicit (specified by equation 19 
below), and other cases without overmodulation are 
discussed. It appears to be advantageous in these 
cases to keep the modulation depth fixed and vary the 
relative amplitude of the two sidebands. The com- 
plications introduced by overmodulation and a way to 
overcome them are discussed in Section 4, 

Only cases in which the upper and lower side- 
bands are initially in phase with the carrier are con- 
sidered, since these are the cases of greatest prac- 



tical importance, and many cases of arbitrary initial 
sideband phase differences are very similar to related 
cases without such phase differences. Theoretical 
results obtained are discussed in Section 5 and 
summarized in Table 1. 

In Section 6 some theoretical results obtained 
are expressed in terms of a readily measurable quan- 
tity and compared with experimental results. Section 
7 gives conclusions. 



2. THE DERIVATION OF S(t) IN THE GENERAL 
TWO-TONE CASE 

An expression for a carrier-wave modulated with a 
signal of n tones is 



E{t) 



cos C0[/ + bi cos (Oit 



i= 1 \ 

\ lai sin (jL>it - bi sin coit j 



) 



cos coot 

(1) 



sin ojot 



(1) 



where 



coo is the (angular) carrier frequency 

(Oi (i = 1, 2„,n) are the (angular) frequencies of the 
modulating tones, 

ai is the amplitude of the lower sideband compon- 
ent at frequency coi 

bi is the amplitude of the upper sideband compon- 
ent at frequency w^ 



and it is assumed ttiat initially both upper and lower 
sidebands are in phase with the carrier* 

The output S{t) to be expected from an ideal en- 
velope detector Is given by 



S{t) 



1 + 



I = 1 



a( cos co\t + hi cos coit 



M^^^w 



i= 1 



ai sin (Lilt ~ hi sin o)it 




(2) 



As already indicated, we shall be concerned only with 
the case in which n = 2. We shall also devote our 
main attention to cases in which the conditions (i) 
(L>± and 0)2 are large compared to their difference, 
(ii) ai = 02 and (ill) h^ = h^ are all fulfilled. When 
n = 2, the modulation depth M is given as a percentage 
by 

M = 100 (ai +02 + ^1 + h>2) (3) 

Consider first the case in which 6Ji = P12, (i^s = 
QI], where P and Q are integers without common 
factor. Then S(t) given by equation (2) is a periodic 
function of t having period 2n/Q., and this is expres- 
sible as a Fourier series 

S(t) = V2 Qo + aicos Q^t + ascos 212^ + . . . (4) 

which has only cosine terms. Suppose that P = 4 and 
Q = 5 (a particular case which is relevant to the com- 
parison of experimental results with theoretical 
results summarized in Table 1 (Section 5) and discus- 
sed in Section 6). A distortion product will have a 
frequency of the form (pwi + gw?), where p and q are 
integers; for the particular case, the distortion prod- 
ucts will therefore have frequencies sH where 



4p + 5g = s 



(5) 



p and q may be positive or negative subject to the 
condition s>0. The. important terms of the particular 
Fourier series for S(t) defined by equations (2) and (4) 
are in general those for which |p| + \q\ is small, 
rather than those for which s is small. To emphasize 
this, it is therefore desirable to specify the coefficient 
of cos s^t in (4) as a^^q rather than as a^. The first 
eleven terms of the Fourier series for S{t) when P = 4, 
Q = 5 are thus 

S(0=V2 cxo,o +cx-i,i cosOi +a.2,2C0s2ni + a2,^iC0s3(li 

+ ai,oCos4nf + OojiCOsSO/ + a-i^2Cos6fli 

+ cx3^»iCos7ni + a2,oCOs8lli + ai,iCos9l]? 

+ ab,2COs10fl/ + (6) 



We expect the most important terms here to be the d.c. 
term and those in 40^, 50i, Q.t, 8(1^, 91)^ and 1011^ 
Similar considerations apply for any other values of 
P and Q. A few specific cases investigated numerical- 
ly indicated that although a smalt change in the value 
of 0)2/^1 may be associated with a large change in 
P, Q and s for a given value of p and q, the change 
in Qp^g is small for fixed oi, a^, fei and ^2. For this 
reason, the case in which 0)2/0)1 is an 'incommensur- 
able' quantity like \/2 which theoretically cannot be 
expressed as the ratio of two integers has been ig- 
nored altogether, and the main interest is in cases 
where P and Q are small; the simplest case is when 
Q~P = 1 so that equals 0J2 ~ wi instead of a frac- 
tion thereof. Now equation (6) is tied to the partic- 
ular values P = 4 and Q = 5, but in the general case 
we can still rearrange (6) so that the terms likely to 
be most important come first, and write 

S(t) = V2 cxo,o + ai^oCOsPfii + ao,iCOsQl]i 

+ a,i,iC0s(Q-P)12i + a2,oCOs2Pni + ai,iCOS(P + Q)n^ 

+ cxo,2COs2Qni + ... (7) 

An obvious way of trying to determine the quantities 
ap^q appears to be to square S{t) in (7) and compare 
this with the squared right-hand-side of (2), equating 
coefficients of cos {pcoi + qco^) t for various values of 
the integers p and q, taken in ascending order of 
\p\ + \q\. The terms actually written down in equation 
(7) are those for which this quantity does not exceed 
2. Terms for which \p\ + \q\ exceeds 2 will of course 
contribute to the coefficients associated with terms 
for which |p| + |g'|<2, but in general we expect such 
contributions to be small, and we start from a case in 
which they are zero in Section 3. 

Since the equations determining the quantities 
Op q are obtained by squaring S{t) in equation (7), the 
equations would equally well be satisfied if every 
^p,q was replaced by -cxp^g* Now in equation (2), 
S{t) is essentially positive, and it is also continuous 
with respect to L But if, in equation (7), ao,o is 
taken as positive, it does not follow that S{t) as given 
by equation (7) will be positive for all values of L 
When the modulation percentage M does not exceed 
100, there is no danger that S{t) given by equation (7) 
can change sign, for in equation (2), S{t) can be re- 
garded as the sum of a fixed vector of magnitude one 
and four rotating vectors of magnitudes ai, a?, /^i and 
fe2- If M<100, then not even if all the four vectors 
are simultaneously in the opposite direction to the 
fixed vector are they capable of cancelling it com- 
pletely. When M>100, such cancellation may or may 
not occur. 

As the comparison of equations (2) and (7) is 
complicated in the general case, we shall apply our 
main analysis to the simpler case where 



ai = 02 = a! 



^2 = a/k 



(8) 



i.e. the two upper sidebands are equai, and tlie two 
tower sidebands are also equal. This case is con- 
sidered in detail in Section 3 when there is no over- 
modulation, and in Section 4 when there is overmodula- 
tion. 



3. THE SPECIAL CASE DEFINED BY EQUATION (8) 

If in the special case defined by equation (8) 
there is no overmodulation, and coefficients ap ^ are 
ignored for \p\ + \q\>2, the Fourier series (7) correctly 
has only cosine terms, and the seven equations giving 
00,0, ai,oi ao,i, oc2,o, ai,i, ^0,2 and a, 1,1 obtained 
by comparing the squares of the right hand sides of 
(2) and (7), are;- 

From the d.c. term - 

1 + 2a^ [1 + (1//2^)] = y4ao,o^ + V2(ai,o^+ ao,i^) 

+ V2(a2,o^ + oi,!^ + ao,2^ + o:-i,i^) (9) 

From the term in oy^it - 

2a [1 +(1//j)] = 

00,001,0 + 0i,o02,0 + 00,101,1 + Oo,iO-i,i (10) 

From the term in co^t ~ 
2a [1 +(1//e)] = 

Oo,oOo,l + Oi^oOi,! + 01,00-1,1+ Oto,l°^, 2 (^1) 

From the term in 2co^t — 

2a^/k -Ob/o02,o +V2ai,o^ + Oi,ia_i^i (12) 
From the term in 2aj^t — 

2a^/k = ob,oab,2 + V2ab,i^ + oi,ia-i,i (13) 
From the term In (a)i + (i>2)? — 
Aa^/k = ob,oOi,i + 01,000, 1 + 02,00.1,1 + 00,20-1,1 

(14) 
From the term in (0)2 -(oi)t — 
2a^[l +(1/fe2)]^ 

ob,oO-i,i +01,000,1 +O2,oai,i +01,100,2 (15) 
Subtracting (13) from (12) - 
oo,o (02,0 -00,2) + V2(ai,o + cxo,i) (ai,o - ao,i) - 

(16) 
Subtracting (11) from (10) - 



(Oi,0O2,0 - 00,100,2) + 

+ (ao,o~a«i,i- ai,i) (ai,o- Oo,i) = (17) 

Equations (16) and (17) are both completely satisfied 
if both 



oi,o = oo,i and 02,0 = oo,2 



(18) 



and we shall assume in what follows that this is the 
only relevant solution, because of the essential sym- 
metry of equations (8). 

If fe = 1, and a<0.25 to prevent overmodulation, 
equations (9) to (15) are exactly satisfied by 

oo,o = 2; oi,o-2a; 02,0 = oi,i = a.1,1 = (19) 

as expected for the nondistortion condition with sym- 
metrical sidebands. (They are also satisfied by 
ao,o = -2 and ai,o = -2a with other ai,j zero as al- 
ready mentioned, but this alternative will not be 
further considered in this Section.) Since S(t) can in 
general be expected to vary continuously when any of 
the quantities involved is slightly altered, we can 
therefore expect ab,o and ai,o to be near these values 
and O2,oi oi,i and a*i,i to be small in neighbouring 
cases. 

If we keep the modulation percentage M constant 
and vary k when there is no overmodulation, a very 
useful initial approximation appears to be 

1 

ab,oJ==J2; ai,oJ==ia(1 +-r) =M/200 (20) 

k 

Substituting this value of ai,o (and ao,i) into (9) and 
neglecting 02,0^, oi,i^, 00,2^ and o_i,i^ gives an 
improved value of 00,0! substituting this into (12) and 
(14) and (15) gives the approximate values of 02,0, 
ai,i arid 0.1,1. Finally substitution into (10) gives 
an improved value of ai,o with which the whole pro- 
cess can be repeated. We have to remember, partic- 
ularly at modulation percentages near 100, that apprec- 
iable contributions may be coming in from higher- 
order terms of frequency pcoi and qco^ where \p\ + 
|q|>2 and that we therefore cannot expect to achieve 
high accuracy in all cases in this way, but the exact 
relation (19) indicates that there are many cases in 
which high accuracy will be achieved. 



4. THE EFFECT OF OVERMODULATION IN SYM- 
METRICAL CASES 

We have already noted that in the case when k was 
1 and a<0.25 (so that there is no overmodulation) the 
exact solution (19) was associated with an ambiguity 
in sign which was fortunately then irrelevant. Phys- 
ically, this kind of ambiguity in sign is associated 
with the possibility of overmodulation, which we must 
now consider. 



\i k = ^ in equation (8), and a<0-25, tlien tine co- 
efficients of S(t) are given by equations (19), so that 

S(t) = 1 + 2a cos co^t + 2a cos co^t (21) 

and tills equation is true wiiatever tiie ratio coi/co2 may 
be. 

If fe = 1 and a>0.25, iiowever, equation (21) lias 
to be replaced by 

S(t) =11 + 2a cos (oit + 2a cos a)2^| (22) 

and equation (19) is no longer valid, because there 
are some ranges of values of t for which 



An Elliott autocode programme has been formulated 
for the evaluation of equations (25) in specific cases. 
This programme has been used to determine the co- 
efficients given in Table 1 (Section 5) for cases with 
and without overmodulation. Most of these calcula- 
tions were made for 0)^/0)2 = 4/5 so that conditions 
(5) applied for determining the important terms for 
which \p\ + |g| is small. Although the calculation has 
to be done in terms of the value of s, the importance 
of the coefficient calculated is primarily controlled by 
the values of p and g. The coefficients given in 
Table 1 (Section 5) embodying results actually cal- 
culated are therefore specified in terms of p and q by 
analogy with equation (7), and not in terms of s. 



(1 + 2a cos co^t + 2a cos o)2t) 

is positive and some for which this expression is 
negative. These ranges are separated by values of 
t for which S(t) is zero, and zeros of S(t) will depend 
upon the value of the ratio ajs/coi as well as upon the 
value of a. 

More generally, when there is overmodulation, 
there is also the possibility that S{t) given by (2) can 
have zeros which are singularities. This means that 
equations analogous to the system (9) - (15) are not 
valid, and that any particular coefficient a^^g will 
depend upon the ratio 0)2/0)1 as well as upon ai, 02, 
foi and &2- There is, however, a well known * sum- 
mation technique'* for obtaining any required Fourier 
coefficient in such cases, when the ratio o)i/ci)2 is 
the ratio of two integers P/Q which have no common 
factor, so that S(t) is periodic with period T, say. 
Suppose that 

S(t) = Uo when ^ = or ^ = T 

(23) 
S(t) = Ur when t = rT/N (r = 1 , 2 



n-d) 



and that we seek an approximation of the form 



R r 



S = 1 



27TSt ^ 27TSt 

as COS — ^ Ps Sin — 



(24) 






where fi<V2(N"1). We shall assume that N is suf- 
ficiently large which means that a^ and ^^ are neg- 
ligible for s>/V-R. Then we can take 

iV-l 

277rs 
cos 

N 

(25) 

JV-1 
2\-^ 2ms 
^s =- 7 l^T sin 

* This is discussed by Clenshaw 'n terms of Chebyshev 
polynomials. 



5. THEORETICAL RESULTS 

Numerical values of the principal coefficients In 
the harmonic analysis of S(t) given by equation (2) 
are tabulated in Table 1 below. It Is assumed that 
0)2/0)1 - 5/4 throughout; only cosine terms are pres- 
ent, k has the same meaning as in equation (8), and 
has one of the values 1. 2, 3, 9 or 00 (the single side- 
band case), /e -= 1 is the equal sideband case, for 
which (21) is an exact solution when there is no over- 
modulation. The principal coefficients are deemed to 
be those listed in equation (7), where for these partic- 
ular calculations means [0)2-0^1)'^ 
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Fig. 1 « Comparison of theoretical and experimental 
values of (0)2-0)1) distortion component 



TABLE 1 

Intermodulation Terms for Two Tones 

ai - 02 = a; ^1=^2 = a/k] M = 200a{k + ^)/k (see equations (2) and (3)). Vaob^o isthe 
d.c. term and (Xp^q is the coefficient of cos (pco^ + q(02)t with o}± = 40, and (02 ^ 511 

(see equation (7)) 



Modulation 
depth (M%) 


k 


a 


0^,0 


«1,0 


ao,i 


«2,0 


«i,i 


CX0,2 


a-1,1 


50 


CX) 


0.25 


2.0633 


0.2440 


0-2440 


-0.0146 


-0.0297 


-0-0146 


0.0318 




9 


0.225 


2.0410 


0-2460 


0-2460 


-0-0096 


-0-0193 


-0.0096 


0.0207 




3 


0.1875 


2.0163 


0.2484 


0-2484 


-0-0038 


^-0077 


-0.0038 


0.0083 




2 


0-1667 


2.0073 


0.2493 


0.2493 


-0.0017 


-0-0034 


-0-0017 


0-0037 




1 


0.125 


2 


0.25 


0.25 














70 


CX) 


0.35 


2.1257 


0.3330 


0.3330 


-0-0265 


-0.0552 


-0-0266 


0-0634 




9 


0.315 


2.0827 


0-3384 


0.3384 


-0.0177 


-0-0364 


-0.0176 


0.0421 




3 


0.2625 


2.0334 


0.3451 


0.3451 


-0-0072 


-0.0147 


-0-0072 


0-0172 




2 


0-2333 


2.0150 


0.3478 


0.3478 


-0-0033 


-0.0066 


-0-0033 


0-0078 




1 


0.175 


2 


0.35 


0.35 














90 


CX) 


0.45 


2.2121 


0.4120 


0.4120 


-0.0387 


^-0836 


-0.0387 


0-1077 




9 


0.405 


2.1427 


0.4228 


0.4228 


-0.0260 


-0.0556 


-0-0260 


0.0737 




3 


0.3375 


2.0597 


0.4376 


0.4376 


-0.0108 


-0.0227 


-0.0108 


0-0316 




2 


0.3000 


2.0274 


0.4441 


0.4441 


-0.0049 


-0.0103 


-0.0049 


0-0146 




1 


0.225 


2 


0.45 


0.45 














100 


00 


0.5 


2.2659 


0-4457 


0.4457 


-0-0431 


-0.0965 


-0.0432 


0.1360 




9 


0.45 


2.1818 


0.4599 


0.4599 


-0.0287 


-0.0639 


-0.0288 


0.0950 




3 


0.375 


2.0785 


0.4804 


0-4804 


-0.0116 


-0-0257 


-0.0115 


0-0425 




2 


0.3333 


2.0368 


0.4902 


0.4902 


-0.0052 


-0.0114 


-0.0050 


0.0203 




1 


0-25 


2 


0.5 


0.5 














110 


00 


0.55 


2.3295 


0.4729 


0.4729 


-0-0443 


-0.1057 


-0-0443 


0-1701 




9 


0.495 


2.2306 


0.4899 


0.4899 


-0.0278 


-0.0678 


-0-0279 


0-1227 




3 


0.4125 


2.1065 


0.5159 


0.5159 


-0.0083 


-0.0235 


-0.0081 


0.0600 




2 


0.3667 


2.0543 


0.5297 


0.5296 


-0-0014 


-0.0072 


-0-0009 


0.0319 




1 


0.275 


2.0033 


0.5470 


0.5469 


0-0024 


0-0026 


0-0026 


0-0025 


130 


00 


0.65 


2.4820 


0-5103 


0.5104 


-0.0377 


-0.1135 


-0.0380 


0.2531 




9 


0.585 


2.3554 


0-5305 


0.5307 


-0.0152 


-0.0629 


-0-0160 


0.1952 




3 


0.4875 


2.1935 


0.5623 


0.5626 


0.0132 


-0.0006 


0.0126 


0.1169 




2 


0.4333 


2.1233 


0.5798 


0.5799 


0.0245 


0-0241 


0.0246 


0.0804 




1 


0.325 


2.0513 


0.6038 


0.6034 


0-0328 


0.0438 


0.0345 


0-0395 


150 


CX) 


0.75 


2.6538 


0.5377 


0.5376 


-0.0288 


-0.1184 


-0.0290 


0-3443 




9 


0-675 


2.4993 


0-5603 


0-5605 


0.0003 


-0.0547 


-0.0013 


0-2772 




3 


0.5625 


2.2996 


0-5960 


0.5970 


0.0390 


0.0266 


0.0361 


0-1856 




2 


0.5000 


2.2117 


0.6160 


0.6170 


0-0558 


0-0610 


0.0533 


0-1423 




1 


0.375 


2.1198 


0-6441 


0.6442 


0.0709 


0.0927 


0.0709 


0.0922 



The coefficients were obtained by means of the com- 
puter programme (based on the first of equations (25)) 
aiready mentioned. The particular ratio coi/a)2 = 4/5 
was chosen, because it was known that experimental 
results were available for the case o)^/2n = 500 Hz 
and 0)2/277 = 621 Hz. The good agreement shown in 
Fig. 1 between experimental and theoretical results 
(and discussed in Section 6) is a useful and encourag- 
ing confirmation of the idea that aparticularcoeff icient 



ap^g is not very sensitive to small changes in the ratio 
0)1 /co2 when p, q, ai, a^, bi and ^2 remain constant, 
when there is no overmodulation. Equations (25) per- 
mit thecalculation of any desired coefficient by means 
of the computer programme, whether there is over- 
modulation or not* But the machine time required for 
such a calculation rises sharply with N, the number of 
values of t for which S{t) must actually be calculated. 
If ccpg is not very sensitive to changes in the ratio 



6 



ws/ct^i, it is clearly advantageous to take tliis ratio 
as (A + 1)/A wliere A is a small integer, if tiiis is at 
all possible, because N can then be relatively small. 

In Table 1 we have confined our attention to 
cases satisfying equation (8), so that there appears to 
be a significant loss of generality. The cases tabu- 
lated, however, appear to be the most important ones 
in practice, and many of the cases which do not 
satisfy equation (8) can be satisfactorily related to 
cases which do satisfy this equation. There is no 
restriction on the validity of equations (25) when 
S(t) is of the form (2), so that any required case can 
be calculated. The above considerations merely 
Indicate that there is no point in actually carrying out 
such calculations in the absence of a specific prac- 
tical need for them. 



6. EXPRESSING RESULTS OBTAINED BY CALCU- 
LATION IN TERMS OF MEASURABLE QUAN- 
TITIES 

A quantity which can be readily measured is 

average modulus of {co2-ci}i) distortion component 



X= 100 



average modulus of total output signal 



(28) 



We require a theoretically determined quantity which 
is comparable with X for cases in which th^re Is no 
overmodulation. 

In all the cases pf interest, the coefficients 
oci,o and oo,i (in the notation of Table 1) are equal, 
and they are much greater than any other coefficient. 
Also, (coq-oji) is a much lower (angular) frequency 
than any other mentioned in connection with Table 1, 
so we have regarded the average modulus of the dis- 
tortion component as approximately equal to the 
average value of |a-i,isin (ci)2-coi)^l over a period 
277/(co2~coi), namely(2/7r)|a-i,i|. The average modulus 
of total output signal is correspondingly taken as 

average value of 

laip cos co±t + ob^i cos co^t + a-i ,i cos (a)2-coi)t| 

(27) 

over the same period. This average modulus of total 
output signal was calculated by means of a computer 



programme when coi = 0-8 q}^, so that the condition for 
the total output signal to be zero reduced to an equa- 
tion of degree 5 in cos (aj2-CL>i)^ which could be 
solved by means of a standard library subroutine. 

Fig. 1 shows experimental values of X for a 
particular receiver (measured at the output of the 
envelope detector) under single sideband conditions, 
for modulation depths between 60% and 100% with 
coi/27r = 500 Hz and ojs/Stt = 621 Hz (broken curve); 
the corresponding calculated values are plotted as a 
continuous curve. The agreement is remarkably good. 
Some further comparisons between theory and experi- 
ment are given in a separate report.^ 



7. CONCLUSIONS 

A complete survey of the general two-tone case is 
perfectly feasible. A computer programme is avail- 
able for analysing any particular numerical case, with 
or without overmodulation. The most important cases 
from the practical point of view — which alone are 
considered in numerical detail here - are those in 
which the two upper-sideband amplitudes are equal, 
the two lower-sideband amplitudes are equal, there 
are no initial differences of phase between corres- 
ponding upper and lower sidebands and the ratio 
ci}'z/(co'z-oj^) is a reasonably small integer. Other 
cases can usually be related closely to the above- 
mentioned special but important cases. The principal 
numerical results obtained are given in Table 1. 

A close agreement between theoretical and meas- 
ured results has been obtained for the quantity X 
specified in equation (26), as indicated in Fig. 1. 
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